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Abstract 

The orchestration of histone modifiers is required to establish the epigenomic status that regulates gene expression during 
development. Whsci (Wolf-Hirschhorn Syndrome candidate 1), a histone H3 lysine 36 {H3K36) trimethyltransferase, is one of 
the major genes associated with Wolf-Hirshhorn syndrome, which is characterized by sl<eletal abnormalities. However, the 
role of Whsci in skeletal development remains unclear. Here, we show that Whsci regulates gene expression through Runt- 
related transcription factor (Runx) 2, a transcription factor central to bone development, and p300, a histone 
acetyltransferase, to promote bone differentiation. Whscl~^~ embryos exhibited defects in ossification in the occipital 
bone and sternum. Whsci knockdown in pre-osteoblast cells perturbed histone modification patterns in bone-related 
genes and led to defects in bone differentiation. Whsci increased the association of p300 with Runx2, activating the bone- 
related genes Osteopontin {Opn) and Collagen type la {Collal), and Whsci suppressed the overactivation of these genes via 
H3K36 trimethylation. Our results suggest that Whsci fine-tunes the expression of bone-related genes by acting as a 
modulator in balancing H3K36 trimethylation and histone acetylation. Our results provide novel insight into the 
mechanisms by which this histone methyltransferase regulates gene expression. 



•0-PLOS I o-^E 



Citation: Lee YF, Nimura K, Lo WN, Saga K, Kaneda Y (2014) Histone H3 Lysine 36 Methyltransferase Whsci Promotes the Association of Runx2 and p300 in the 
Activation of Bone-Related Genes. PLoS ONE 9(9): el06661. doi:10.1371/journal.pone.0106661 

Editor: Mary Bryk, Texas A&M University, United States of America 

Received March 20, 2014; Accepted August 6, 2014; Published September 4, 2014 

Copyright: © 201 4 Lee et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its 
Supporting Information files. 

Funding: Funding came from MEXT KAKENHI Grant Nos. 23710224 and 24116514 and the Osaka University Program for the Support of Networking among 
Present and Future Researchers. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 

* Email: nimura@gts.med.osaka-u.ac.jp 



Introduction 

Histone modifications are crucial in ensuring proper embryonic 
development [1] [2]. For example, Wolf-Hirschhorn syndrome 
(WHS) is a syndrome closely associated with epigenetics. WHS is 
defined by craniofacial dysgenesis, growth delay, mental retarda- 
tion, and heart malformations, among other characteristics [3] [4] , 
and is caused by sub-telomeric deletions on the short arm of 
chromosome 4pl6.3. Whsci (WSH candidate 1, also known as 
NSD2 or MMSET), the gene that encodes the histone 3 lysine 36 
(H3K36) trimethyltransferase Whsci [5], is deleted in every case of 
WHS [3], and this deletion is necessary for the occurrence of 
WHS [6]. Because WHS patients show skeletal abnormalities, 
such as sternal hypo-ossification [7], delayed growth, and 
craniofacial abnormalities [3] [8], it has been concluded that 
Whsci is involved in bone development. However, its function in 
skeletal development remains unknown. 

Runt-related transcription factor (Runx) 2 plays a critical role 
during the differentiation of mesenchymal cells into osteoblasts [9] . 
Runx2 belongs to the Runx family of genes, which contain the 
conserved Runt DNA binding domain and are key regulators 
during development. Gene knockout models have helped elucidate 
the key functions of the Runx genes; Runxl is associated with 
hematopoiesis, Runx2 with bone formation, and Runx3 with 
cytotoxic T-ceU development [10]. Runx genes have many other 



roles in addition to their major functions and can be co-expressed 
in the same tissue. Although Runxl plays a key role in 
hematopoiesis, it also cooperates with Runx2 [1 1] [12] in sternum 
development and in the early stages of mesenchymal differenti- 
ation [13]. Runx genes function as DNA binding transcription 
factors and interact with co-repressors or co-activators in the 
suppression or activation of genes. Runx2 represses the osteocalcin 
gene [14] when it interacts with the co-repressor histone 
deacetylase (HDAC) 3. In contrast, Runx2 activates matrix 
metalloproteinase-13 gene [15] when bound to p300, a co- 
activator histone acetyltransferase (HAT). The co-regulation of 
histone modification-related enzymes and bone-related transcrip- 
tion factors, including Runx2, has been found to be crucial in the 
regulation of bone cell differentiation genes [16]. 

Whsci contains the SET (suppressor of variegation, enhancer of 
zest and trithorax) domain, which specifically catalyzes histone 3 
lysine 36 methylation (H3K36me). H3K36me can exist in non- 
methylated (meO) or mono-, di-, or tri- methylated forms (mel, 
me2, or me3, respectively). With respect to the methylation status 
of H3K36, H3K36me3 has been highly correlated with tran- 
scribed genomic regions from transcription start sites to transcrip- 
tion termination sites [17] [18]. Other HMTs that catalyze 
H3K36me3 include HYPB/SETD2 and NSDl [19] [20] [21]. In 
yeast, SET2, the ortholog of HYPB/SETD2, methylates H3K36 
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in transcribed regions. Througli its interaction witti the elongating 
RNA polymerase II, cryptic transcription, in which transcription 
begins in the middle of the coding region, is inhibited via the 
suppression of histone turnover [22] [23]. HYPB/SETD2 acts in 
conjunction with RNA polymerase II, whereas Whscl interacts 
with transcription factors to prevent the overactivation of genes 
during heart and embryonic stem cell development [5]. 

In the present study, we observed that Whscl associates with 
Runx2 in the regulation of bone-related genes, ensuring proper 
skeletal development. In particular, Whscl fine-tunes the expres- 
sion of bone-related genes through H3K,36me3 and histone 
acetylation by promoting the interaction between Runx2 and 
p300. These findings suggest an essential role for Whscl with 
Runx2 and pSOO in the activation of bone-related gene expression 
during bone development. 

Materials and Methods 

Ethics 

We used mice in accordance with protocols approved by the 
Ethics Committee for Animal Experiments of the Osaka 
University Graduate School of Medicine. 

Cell culture and mice 

Murine C3H lOTl/2 mesenchymal ceUs [24] and MC3T3-E1 
murine pre-osteoblasts [25] were used in the in vitro studies. C3H 
lOTl/2 mesenchymal cells were grown in Dulbecco's modified 
Eagle's medium (DMEM) (Nacalai Tesque) supplemented with 
10% fetal bovine serum (FBS), 100 U/ml penicillin, and 0.1 mg/ 
ml streptomycin. MC3T3-E1 cells were cultured in a-modified 
Eagle's medium (a-MEM) (GIBCO Catalogue: 12571-063) 
supplemented with 10% FBS, 100 U/ml penicillin, and 0.1 mg/ 
ml streptomycin. Osteoblast differentiation was induced by the 
addition of 300 ng/ml bone morphogenetic protein (BMP)-2 
(PEPROTECH) to the culture medium. 

Whscl~^~ embryos were obtained from Whscl^^~ mice that 
had been backcrossed to C57BL/6 over 10 times. The targeting 
strategy was conducted as described in [5]. 

Skeletal preparations 

Embryos were fixed in 2()"'o formalin for 2 days at room 
temperature and subsequc-ntly e\ iscerated. After dehydration in an 
increasing ethanol series up to 100% ethanol for 30 minutes each, 
the embryos were bleached with Dent's fixative (5:36:9 ratio of 
aqueous H202:MeOH:DMSO, respectively) overnight. The em- 
bryos were then washed with a decreasing ethanol series, down to 
70% ethanol, for 30 minutes each. Suhsc-quently, tlu; embryos 
were stained with 0.1% Alizarin red at 37°C overnight. The 
embryos were again dehydrated in an ethanol series, up to 100% 
ethanol, for 30 minutes each, followed by staining with Alcian blue 
solution for 2 days at 37°C. The embryos were then washed and 
placed in xylene for 2 hours to remove fats. Then, the embryos 
were cleared with a trypsin solution (30 ml: 70 ml:l g ratio of 
sodium tetraborate saturated solution:distilled water:trypsin, 
respectively) at 37°C and kept in 80% glycerol. 

Micro CT imaging of newborn was obtained by in vivo micro 
X-ray CT system R_mCT2 (Rigaku). 

Alkaline phosphatase assay 

MC3T3-E1 cells plated in 24-well plates were washed twice 
with Tris-bufiered saline (TBS), scraped into 0.3 ml of 0.5% NP- 
40 containing 1 mM MgClj and 10 mM Tris (pH 7.5), and 
homogenized three times for 10 seconds each with a multi-bead 
shocker (Yasui Kikai). Then, the cell lysates were centrifuged for 



5 minutes at 2000 rpm, and the supematants were used for the 
enzyme assay. Alkaline phosphatase activity was assayed using a 
commercial kit (Wako Pure Chemicals Industries) with p- 

nitrophenylphosphate as a substrate. The enzymatic activity was 
expressed in units/|J,L defined by the release of 1 nmol of p- 
nitrophenylphosphate per minute at pH 9.8 and 37°C. This 
enzymatic activity was then normalized to the protein content of 
the sample. The protein content was determined using a 
bicinchoninic acid (BCA) protein assay kit (Thermo Scientific) 
with bovine serum albumin (BSA) as the standard. 

Plasmid construction 

To create Runxl and Runx2 plasmids for the overexpression 
experiments, mouse Runxl and Runx2 were first amplified using 
polymerase chain reaction (PGR) from cDNA templates, which 
were reversed transcribed from the total mRNA of C3H lOTl/2 
cells. Runxl, Runx2, and Runx2 deletion mutants were 
introduced into the pGAGIP-myc and pCAGIP-HA expression 
vectors using Gateway Technology (Invitrogen). The Whscl 
expression constructs have been previously described [5]. For 
the luciferase promoter assays, the Osteoponlin (Opn) promoter 
(—1010 to +146) and the Collagen type la (Co/7a7)-luciferase 
construct (—2320 to -1-183) were generated by PGR using mouse 
genomic DNA, and the resulting fragments were cloned into the 
pGL3-basic vector (Promega). The primer sequences are shown in 
Table SI. 

Transient plasmid transfections and targeted gene 
knockdown 

FuGENE HD (Promega) was used to transfect plasmids into 
GH3 lOTl/2 cells according to the manufacturer's instructions. A 
mixture of 35 jil of FuGENE HD and 10 ng of DNA was used to 
transfect the cells seeded in a 10-cm dish. The Neon (Invitrogen) 
electroporation method was used to transfect plasmids and small 
interfering RNAs (siRNAs) at a final concentration of 100 nM into 
MG3T3-E1 cells. The following conditions were used for Neon: 
pulse voltage, 1600 V; pulse width, 10 ms; and pulse number, 3. 

The siRNA method was used to knockdown the expression of 
targeted genes in MG,3T3-E1 cells. Runxl, Runx2, Whscl, or 
control siRNA (Sigma-Aldrich) were transfected into MG3T3-E1 
cells using tiu- Nc-on Transfection System (Invitrogen) as mentioned 
above. MC3T3-E1 cells were cultured for 1-2 days in a 12-well 
plate at a cell density of 1 xlO^ cells/well. The following siRNAs 
were obtained from MISSION (Sigma-Aldrich): siRunxl.6068 
(SASI_Mm02_00306068), siRunxl.6069 (SASI_Mni02_003060 
69), siRunxl.6070 (SASI_Mm02_00306070), siRunx2.4517 (SA- 
SI_Mm01_00044517), siRunx2.4518 (SASI_Mm_00044518), siW 
hscl.5204 (SASI_Mm02_00295204), and siWhscl.5207 (SA- 
SI_Mm02_00295207). 

Cell fractionation 

To obtain nuclear extracts for analyses, MC3T3-E1 (IxlO' 
ceUs/sample) and C3H lOTl/2 (3.7x10'' cells/sample) cells were 
extracted in 500 |J.l of nuclear isolation buffer (NIB) containing 
0.15% NP-40 and EDTA-free 'complete' protease inhibitor 
cocktail (Roche) as previously described [26] . Cytoplasmic soluble 
proteins were separated from the nuclei by centrifugation at 
15,000 rpm for 5 minutes at 4°C. Nuclear pellets were further 
treated with 1200 U of micrococcal nuclease (TaKaRa) in 250 (xl 
of nuclear isolation buffer containing 200 mM NaCl at 25°C for 
15 minutes. After incubation on ice for 10 minutes, 10 mM 
EDTA was added, and the samples were incubated on ice for 
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10 minutes. The nuclear extracts were separated from the pellets 
by centrifugation at 15,000 rpm for 5 minutes. 

Co-immunoprecipitation 

C3H lOTl/2 cells were used to study the interactions between 
endogenous p300 and the overexpressed Runx and Whscl 
proteins. Nuclear extracts were isolated from C3H lOTl/2 cells 
24 hours after transfection. Nuclear extracts were first pre-cleared 
with rabbit IgG AC (Santa Cruz Biotechnology) and then 
incubated overnight with HA-tagged beads to bind the Runx- 
HA-tagged proteins. The beads were washed 4 times using nuclear 
isolation buffer (NIB) containing 200 mM NaCl. 

RNA extraction and real-time PCR 

To analyze gene expression, total RNA from embryonic stage 
18.5 (E18.5) mouse bone tissues and MC3T3-E1 cells were first 
extracted using Isogen (Nippon Gene). For bone tissues, homog- 
enization was performed before extraction. Tissues were homog- 
enized three times at 2000 rpm for 10 seconds in 1 ml of Isogen 
and then centrifuged at 15,000 rpm for 5 minutes to remove cell 
debris. The supernatant was collected, and RNA extraction was 
performed according to the manufacturer's instructions. In each 
experiment, an equal quantity of total RNA was reversed- 
transcribed into cDNA using Superscript III (Invitrogen). Quan- 
titative PCR was performed using SYBR Premix Ex Taq 
(TaKaRa) and Real-time PCR Master Mix (TOYOBO). Amph- 
fications were performed in a CI 000 thermal cycler (Bio-Rad). 
The quantified mRNA levels were normalized to the expression of 
Rplp2 (bone tissues) and 18S rRNA (MC3T3-E1) mRNA (as an 
internal control). TaqMan probes and primer pairs specific for 
Osleoponlin (Sppl) (TaqMan Mm00436767_ml), Osteocalcin 
{Bglapl) (TaqMan Mm03413826_mH), Alkaline phosphatase 
(Mm00475381_ml), Collagen type la (TaqMan Mm00801666_ 
gl), Pvunxl (TaqMan Mm01213405_ml), and Whscl (TaqMan 
Mm01211 104_ml) were purchased from Applied Biosystems 
(Foster City, CA). The primers for Runx2, CollOal, Rplp2 
(Ribosomal protein, large P2), and 18S {18S rRNA) are shown in 
Table 81. Assays were performed in triplicate within a single 
experiment, and the average values were used for further analysis. 

Pre-osteoblast cell count 

MC3T3-E1 cells were counted at 36 hours after transfection 
using a cell counter (Beckman Coulter Z) or a TC-10 Automated 
Cell Counter (Bio-Rad Laboratories). 

Luciferase reporter assay 

The Opn or Collal luciferase reporter constructs and the 
expression constructs were co-transfected into C3H lOTl/2 cells 
that were seeded at a density of 2x10* cells/well in 24-^\-cll plates 
one day prior to transfection. The ratios of the expression 
constructs were as follows: for the luciferase reporter construct, the 
Renilla plasmid ratio, as an internal control, was 5:4:1, and 0.5 (ig 
of total DNA was transfected per well. The transfection reagent/ 
DNA complex for 10 wells was made from 5 |J,g of DNA and 
17.5 \i\ of FuGENE HD (Roche) in 250 \i\ of OPTI-MEM. A total 
of 25 |J,1 of the resulting complex was added to each well. C646, a 
p300 inhibitor (Sigma-Aldrich), was dissolved in DMSO and was 
added to the appropriated groups 1 2 hours after transfection at a 
concentration of 45 nM. After 24 hours of incubation, luciferase 
activity was measured using the Dual Luciferase Reporter Assay 
system (Promega). 



Western blot analysis 

Trypsinized MC3T3-E1 cells were washed once with PBS and 
dissolved in 2X sample buffer (1 x 10' cells/ 10 |il). To denature the 
proteins, the samples were vortexed for 5 minutes and incubated 
at 100°C for 5 minutes, and this procc-ss was repeated 3 times. 
Immunoprecipitation samples were also analyzed by western blot 
analysis. Protein samples (20 ^1 per lane) were separated by SDS- 
PAGE and transferred onto a Hybond-P PVDF membrane (GE 
Healthcare). The following antibodies were used for western blot 
analysis: anti-histone H3 (1:1000, abl791. Abeam), anti-Runxl 
(1:500, 8529, Cell Signaling Technology), anti-Runx2 (1:1000, 
8486, Cell SignaUng Technology), anti-p300 (1:100, sc-585, Santa 
Cruz Biotechnology), and anti-Whscl (1:200) [5]. AH of the 
primary antibodies were detected using secondary antibodies 
conjugated to horseradish peroxidase (anti-rabbit IgG, GE 
healthcare). A ReliaBLOT (Bethyl Laboratories) kit was used to 
reduce the IgG heavy chains in the IP samples. To strip the blots, 
the membranes were incubated with WB Stripping Solution 
Strong (Nacalai Tesque) for 15 minutes at room temperature. 

Chromatin immunoprecipitation (ChIP) qPCR 

MC3T3-E1 cells were transfected with control siRNA or siRNA 
directed against Whscl and were cultured in BMP-2 for 6 days. 
Nuclei was isolated by NIB containing NP40 from 0.5 xlO' cells, 
then chromatin was digested at 25°C for 30 min by 4.8 U ml ' 
micrococcal nuclease (TaKaRa) in 250 |j.l of NIB contained 
400 mM NaCl as described previously [26]. The nuclear extract 
was diluted to 200 mM NaCl by 250 ^1 of NIB contained 0 mM 
NaCl. Chromatin was immunoprecipitaed with anti-Runx2 
(2.5 8486, CeU Signahng Technology), anti-H3K36me3 
(5 ng, ab9050, abeam), and anti-H4ac (5 |ll, 06-866, MiUipore) 
antibodies, then purified using a QIAquick PCR purification kit 
(Qiagen) [27]. THUNDERBIRD'^' SYBR® qPCR mix (TO 
YOBO) was used for real-time PCR analysis. The primer 
sequences are shown in Table SI. 

Statistical analysis 

Statistical anah'scs wviv pc'rform(;d using JMP 9 software (SAS 
Institute). Data were compared using the two-tailed Student's /-test 
for pairwise comparisons against control and the Tukey-Kramer 
HSD test for multiple comparisons. Significance was defined as 
P<0.05. 

Results 

Skeletal abnormalities in the Whscl knockout mice 

To study the aberrations of the skeletal system, whole skeletal 
structures of E18.5 wild type and Whscl ' embryos were fixed 
and stained with Alizarin red to stain the calcified bones and 
Alcian blue to stain the cartilage. Whscl~^~ embr^'os had 
significant disruptions of ossification in cranial bone elements, 
including the occipital bone and periotic bone, as denoted by 
arrows in Fig. lA. An enlarged top view of the cranial bone 
indicated a severe lack of ossification in the occipital bone 
(Fig. IB). In addition to the abnormal cranial bone ossification, the 
embryo had insufficient ossification in the sternum (Fig. lA and 
IC). We next examined the alkaline phosphatase (Alp) activity as 
an early marker for osteoblast differentiation. Whscl deficiency 
significantiy decreased the Alp activity in the occipital bone and 
sternum (Fig. ID). To examine the mineralization aberrations of 
Whscl ' bone, Whscl ^ newborn was analyzed by micro X-ray 
CT system, because Whscl knockout is lethal after birth [5]. 
Whscl ' newborn showed mineralization deficiency in occipital 
bone, sternum, and clavicles (Fig. IE). These results suggest that 
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Wild Type Whsd-/- Wild Type Whsc1-l- 




Figure 1. Skeletal abnormalities in E18.5 Whscr^^ embryos. (A-C) Calcified tissues from wild type and Whscl~'~ embryos were stained red 
with Alizarin red, and the cartilage was stained blue with Alcian blue. (A) Whole skeletal structure. Arrows represent the areas of the skeletal structure 
lacking ossification. (B) Top view of the cranial bones. Arrows denote the extensive lack of ossification in the occipital bone. (C) Ribs with sternum. (D) 
Alkaline phosphatase assay on bone samples from wild type and Whscl^^^ embryos. *, P<0.05. (E) IVIicro-CT imaging of mineralized bones in 
postnatal 0 day wild type and Whscl^''^ newborn. Arrow heads indicate sternum. Arrows indicate occipital bone or clavicles. 
doi:1 0.1 371/journal.pone.01 06661 .gOOl 
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Whscl is involved in the regulation of skeletal development, 
notably in the sternum and occipital bone. 

Whscl deficiency decreased the expression of bone 
development-related genes 

To understand the genetic basis of the ossification deficiencies 
observed in the occipital bone and sternum of Whscl ^ embryos, 
we examined the expression levels of difiFerent genes in several 
bone tissues of E18.5 Whscl~^~ embryos. Whscl deficiency did 
not change the expression levels of the Runxl and Runx2 genes, 
the critical transcription factors for bone development, or the 
bone-related proteins Opn, Osteocalcin (Ocn), Collal, and Alp in 
the frontal bones (Fig. 2A) and limbs (Fig. 2B). In the occipital 
bone, Whscl deficiency decreased the expression of Opn and 
Collal but not Runxl and Runx2 (Fig. 2C). The sternum of the 
Whscl ^ embryos exhibited a drastic decrease in the expression 
of both transcription factors and bone-related proteins (Fig. 2D 
and Fig. SI). These results suggest that Whscl plays critical roles in 
bone differentiation in the sternum and occipital bone. 

Whscl is involved in osteoblast differentiation but not in 
proliferation 

To study the role of Whscl in bone diHerentiation, we repressed 
Whscl expression in pre-osteoblastic MC3T3-E1 cells. Bone 
differentiation was induced by the administration of BMP-2, which 
is known to induce osteoblast differentiation [28]. Cells cultured 
for 3 days showed a minimal increase in osteoblast difiFerentiation- 
related gene expression compared with cells cultured for 6 days 



(Fig. 3A, 3B, and Fig. S2). Therefore, MC3T3-E1 cells started 
differentiating after 3 days in our method. Significant decreases in 
the levels of the differentiation-related markers Alp (Fig. 3A) and 
Ocn (Fig. 3B) were observed in the Whscl-, Runxl-, and Runx2- 
depleted cells, indicating that these proteins are important for 
osteoblast differentiation. The decrease in gene expression was 
especially significant in the Runx2-knockdown cells, which is 
consistent with the report that Runx2 directiy regulates osteoblast- 
specific genes, such as Opn [29]. Whscl -knockdown cells exhibited 
a substantial difference in Alp activity after 6 days of BMP-2 
treatment (Fig. 3C). Histochemical Alp staining further supports 
the above data, with lower Alp staining in the Whsc 1 -knockdown 
cells compared with the controls (Fig. 3D). Whsc 1 -knockdown also 
decreased Collal, Opn, and Idl expression in MC3T3-E1 cells 
after 6 days of BMP-2 treatment (Fig. 3E). Taken together, the 
significant decrease in osteoblast differentiation-related genes and 
the decrease in Alp production observed in Whsc 1 -knockdown 
cells indicate that Whsc 1 is necessary for osteoblast differentiation 
at least in occipital bone and sternum. 

Next, to exclude the possibility that Whscl affects the 
proliferation of pre-osteoblast cells, we examined whether Whsc 1 
knockdown decreases cell proliferation. Whscl knockdown did not 
decrease the expression of Runxl or Runx2 in pre-osteoblast 
MC3T3-E1 cells, nor did depletion of Runxl or Runx2 affect 
Whscl expression (Fig. 3F). No significant differences in the 
relative cell numbers were observed between the control and the 
Whscl-, Runxl-, or Runx2-depleted cells after 6 days of BMP-2 
treatment (Fig. 3G), indicating that the Whscl and Runx genes do 





Runxl Runx2 Opn Ocn CoMal Alp Whscl 



■Wild Type 
■ Whsd-/- 



Sternum and Sternebrae 



^1 



Runxl Runx2 Opn Ocn CoHal Alp Whsd 



Figure 2. Whsd deficiency decreases the expression of bone tissue-related genes in El 8.5 occipital bones, sternum and sternebrae. 

(A-D) Gene expression in Whscl^^^ El 8.5 bone tissues was measured by qRT-PCR. Averages from six embryos of each genotype group are shown. 
Values were normalized to Rplp2. Error bars represent the s.d. *, P<0.05. **, P<0.001. 
doi:1 0.1 371/journal.pone.01 06661 .g002 
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Figure 3. Whsci deficiency Inhibits osteoblast differentiation but not cell proliferation. (A, B) Alp (A) and Ocn (B) mRNA expression during 
bone differentiation was measured by qRT-PCR. IVIC3T3-E1 cells were transfected with the Indicated sIRNAs and differentiated by Incubation with 
BIVlP-2. Values were normalized to the expression of 18S. (C, D) Alp activity was assessed using an alkaline phosphatase assay (C) or through Alp 
staining (D) In MC3T3-E1 cells transfected with the Indicated sIRNAs and In the presence or absence of BMP-2. (E) mRNA expression during bone 
differentiation was measured by qRT-PCR. MC3T3-E1 cells transfected with the Indicated sIRNAs were differentiated by 3 days (BMP3) and 6 days 
(BMP6) Incubation with BMP-2. Values were normalized to the expression of 18S. (F) MC3T3-E1 cells transfected with the Indicated sIRNAs were 
analyzed by western blot analysis with the Indicated antibodies. (G) Cell proliferation of MC3T3-E1 cells transfected with the Indicated sIRNAs. MC3T3- 
E1 cells were differentiated by 6 days Incubation with BMP-2. The value of wild type was set to 1 .0. For all panels, Error bars represent the s.d. *, P< 
0.05. **, P<0.01. 

dol:1 0.1 371 /journal.pone.01 06661 .g003 



not regulate cell proliferation before difTereiitiation. Therefore, 
these results suggest that Whsc 1 , Runx 1 , and Runx2 are important 
for osteoblast differentiation but not for the proliferation of pre- 
osteoblasts. 

Whsci interacts with the skeletal transcription factor 
Runx2 and its co-activator p300 

Next, we examined whether Whsci was associated with Runx2. 
Whsc 1 co-immunoprecipitated with Runx2 as well as the dN and 



drunt truncated proteins, but not it did not co-immunoprecipitate 
with the drunxl truncated protein (Fig. 4A and 4B). This suggests 
that Runx2 interacts with Whsci through the runxl domain and 
hence Runx2 may associate with Whsci in the promotion of 
osteoblast differentiation. The histone acetyltransferase (HAT) 
p300 is known to interact with Runxl [30] and Runx2 [3 1] during 
transcription activation. In heart development, Whsci led to a 
decrease in the H3K36me3 levels, and increased gene expression 
was obsei-ved [5]. However, the expression of bone-related genes 
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was decreased in Whscl embryos (Fig. 2C and 2D). Thus, 
Whscl could recruit p300 to the Runx complex. Indeed, Whscl 
overexpression increased the association of p300 with Runx2 
(Fig. 4C). Although interactions between p300, Runxl, and 
Whsc 1 were detected, no increase was observed in the association 
between p300 and Runxl (Fig. 4C). Non-specific binding of 
Whscl to control IgG was found, but p300 was not precipitated at 
the same time, suggesting that Whscl that was not involved in 
Runx2 complex was non-specifically precipitated with control IgG 
(Fig. 4C). Therefore, these results suggest that Whscl promotes 
the association of p300 and Runx2 but not the association of p300 
and Runxl. 



Whscl is crucial for Opn and Collal gene activation 

To reveal the function of the association of the Runx proteins 
with Whscl and p300 in promoting osteoblast difiFerentiation, we 
examined whether these factors regulated Opn and Collal 
promoter activity. When Whscl and Runx2 were overexpressed, 
Opn and Collal promoter activity was significantly increased 
(Fig. 4D and 4E). In contrast, overexpression of Runxl and 
Whscl did not increase Opn and Collal promoter activity 
(Fig. 4D and 4E). Furthermore, the p300 inhibitor C646 repressed 
the activation of the Opn and Collal promoters following Whscl 
and Runx2 overexpression (Fig. 4D and 4E). A mutation in the 
SET domain of Whscl (H1143G), which disabled the histone 
methyltransferase (HMT) ability of Whsc 1 (Fig. 4F) [5] , was also 
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used to investigate the effect of H3K36me3 loss on promoter 
activation. Overexpression of H1143G and Runx2 caused 
stronger promoter activity compared with wild type Whscl and 
Runx2 overexpression (Fig. 4D and 4E) in a way that is consistent 
with the report that the H3K36me3 activity of Whscl represses 
overactivation of Pdgfra gene in heart development [5] . These 
data suggest that Whscl, together with Runx2 and p300, is 
important for promoter activation, and its HMT activity 
modulates transcription level of the Opn and Collal genes. 

Whscl deficiency decreases Runx2 binding and histone 

H4 acetylation in target genes 

To confirm the effect of Whscl on the Opn and Collal 
genomic regions, Runx2 and the histone status at the Opn and 
Collal loci were examined using ChIP qPCR. First, We 
examined Runx2 binding to these genes in wild type MC3T3- 
El cells in 6 days culture with BMP2. Runx2 bound to Opn and 
Collal loci, compared with Ocn promoter including Runx2 
binding sites [32] as positive control and Hlfoo, an oocyte specific 
gene, as negative control (Fig. 5A). H3K36me3 was highly 
enriched at Collal exonl region and moderately enriched at 
other regions compared to Hlfoo (Fig. 5A). The acetylation of H4 
(H4ac) was enriched at Opn promoter, Collal promoter, and 
Collal cxonl regions (Fig. ,^A). Whscl knockdown significantly 
decreased Runx2 binding to the promoter and first exon regions of 
Opn, the first exon oi Collal, and the promoter of Ocn genes in 6 
days BMP-2-treated MC3T3-E1 cells (Fig. 5B). H3K36me3 levels 
in the first exon regions were also decreased in Whscl -depleted 
MC3T3 cells (Fig. bC and 5D). H4ac was measured to assess the 
effect of p300. H4Ac levels were decreased at the promoter oiOpn 
and Collal in Whscl -depleted cells (Fig. 5E and 5F). These 
results suggest that Whscl is important for recruiting Runx2 and 
p300 to the promoters of bone-related genes. 

Discussion 

Our results provide evidence for the role of Whscl in skeletal 
development. We demonstrated that Whscl, a H3K36 trimethyl- 
transferase, is crucial in modulating the transcriptional activation 
of bone-related genes, such as Opn and Collal, through its 
associations with Runx2 and p300. In El 8.5 Whscl ^ embryos 
and newborns, significant disruptions in ossification were observed 
in the sternum (Fig. lA, IC, and IE) and some cranial bone 
elements, particularly in the occipital bone (Fig. lA, IB, and IE), 
but not in frontal bone and limb (Fig. lA and IE). Setd2 is also 
reported to have the activity of H3K36me3 [20] [33]. Because 
Whscl and Setd2 have the same enzymatic artivity, they might 
have redundant function in the regulation of gene expression 
during bone development. Ossification disruption by Whscl 
deficiency in selective bones might be due to Setd2 rescuing 
Whscl deficiency in frontal bone and limb. The cranial bone 
defects observed in the Whscl ^~ embryos may explain the 
craniofatdal defects observed in WHS patients. The phenotype of 
Whscl embryos affirms a link between Whscl and skeletal 
development. Whsc 1 deficiency significantiy decreased Alp activity 
(Fig. ID) and Opn and Collal gene expression (Fig. 2C and 2D) 
in the sternum and occipital bone. Furthermore, Whscl knock- 
down in MC3T3-E1 cells also disrupted the progress of bone 
differentiation without suppressing cell proliferation (Fig. 3). 
However the gene expression of Ocn, a marker of mineralization, 
was not significantiy changed in Whscl~^~ occipital bone and 
sternum (Fig. 2C and 2D), suggesting that Whscl deficiency may 
also affect early commitment or the numbers of osteoprogenitors. 



These results indicate that Whscl is essential to ensure proper 
ossification through the regulation of these bone-related genes. 

Skeletal development is an intricate process that requires 
stringent control of gene activation and suppression by a variety 
of transcription factors [34]. Because Whscl interacts with heart 
transcription factors to ensure proper heart development [5], 
Whscl most likely associates with bone transcription factors. 
Therefore, because Runx transcription factors play pivotal roles in 
skeletal development, we focused on the association of Whscl with 
the Runx genes. Because the mRNA expression of Runxl and 
Runx2 in the sternums of Whscl~'~ embryos (Fig. ID) was 
decreased, we sought to determine whether Whscl directiy 
regulated the Runx genes. Whscl knockdown did not substantially 
alter Runx protein expression (Fig. 3F). This result suggests the 
possibility that defects in bone differentiation caused by Whscl 
deficiency are not mediated by a decrease in Runx protein 
expression. Therefore, Whscl may instead act as a co-regulator of 
Runx proteins in a manner similar to the association between 
Whscl and Nkx2-5 observed during heart development [5]. 

Because Whscl may be associated with Runx in the promotion 
of osteoblast differentiation, we sought to confirm that this 
complex is involved in the transcriptional control of bone-related 
genes. We initially thought that Runxl may be more closely 
associated with Whscl because conditional Runxl knockout 
embryos [1 1] exhibited a similar lack of ossification in the sternum 
and occipital bone. However, our results indicated that although 
Runxl was associated with Whscl (Fig. 4C), the combination of 
Runxl with Whscl did not significantly enhance Opn and Collal 
promoter activity, unlike the significant enhancement observed for 
Runx2 and Whscl (Fig. 4D and 4E). Thus, our results indicate 
that Whscl regulates bone-related gene expression in conjunction 
with Runx2. 

Although H3K36me3 is a mark of transcribed genomic regions 
and has the function to repress overactivation of transcription in 
active genes [5], we showed that Whscl is a positive regulator of 
bone-related genes (Fig. 4D and 4E). These data imply that Whscl 
requires co-activators for the activation of promoters. Histone 
acetyltransferases such as p300, are involved in transcription 
activation [35]. In addition, p3()() recruitment to Runx2 has also 
been reported to play a role in the transcriptional activation of 
bone-related genes [31], suggesting that HATs, such as p300, may 
be associated with Whscl and Runx2. Our co-immunoprecipita- 
tion results indicated an association between Whscl, Runx2, and 
p3()() (Fig. 4C). Moreover, Whscl promoted the association 
between Runx2 and p3()() (Fig. 4C). The requirement of p300 
for the activation oiOpn and Collal was further demonstrated by 
the decreased promoter activity after the addition of the p300 
inhibitor C646 (Fig. 4D and 4E). In addition, the level of H4Ac, 
which is deposited by p300, was correlated with Whscl (Fig. 5D 
and 5E). Although H3K36me3 catalyzed by the HMT activity of 
Whscl was not required for recruiting Runx2 and p300 to 
promoter and for activating the bone differentiation genes, its loss 
led to the overactivation of the promoters (Fig. 4D and 4E). Thus, 
H3K36me3 catalyzed by Whscl allows for repressive mechanisms 
to prevent overactivation of target genes. Therefore, we propose 
that the Opn and Collal promoters are transcriptionally regulated 
during osteoblast differentiation. Whscl acts as a scaffold protein 
to recruit p300, thus promoting the association between Runx2 
and p300. The activation of the Opn and Collal genes occurs in 
the presence of this activation complex, with H4Ac catalyzed by 
p300. Whscl also deposits H3K36me3 through its HMT activity 
on the first exon region (Fig. 5C and 5D) and this represses 
overactivation. Together with the activating acetylation marks of 
p300, the proper level of promoter activation occurs. Runx2 was 
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reported to regulate bone difTerentiatioii in dose-dependent 
manner, suggesting that a critical threshold concentration of 
Runx2 is required for adequate transcription level [36] [37]. Thus, 
the function of Runx2 might be finely regulated for precise gene 



expression. Runx2-associating factors including Whsci and p300 
might be critical for Runx2 activity to activate bone differentiation 
genes, because Whsci knockdown decrease not only Runx2 
binding to the target regions of Opn and Collal genes but also 



PLOS ONE I www.plosone.org 



9 



September 2014 | Volume 9 | Issue 9 | el 06661 



Whsci Fine-Tunes Bone-Related Genes 



histone acetylation on the promoter of these genes (Fig. 5). 
Furthermore, the embryos with 50% of the wild type Runx2 level 
show severe defect in occipital bone, compared with frontal bone 
[36], which is similar to Whscl ^ embryos. These suggest that 
Whscl fine-tunes the transcription of the Opn and Collal genes 
with Runx2 and p300. 

In conclusion, we propose a possible cross-talk pathway between 
Whscl, Runx2, and p300 in the transcriptional activation of bone- 
related genes, leading to proper skeletal development at least in 
occipital bone and sternum. 

Supporting Information 

Figure SI GollOalexpression in Whscl~^~ E18.5 bone 
tissues was measured by qRT-PCR. Averages from six 

embr)'Os of each genotype group are shown. Values were 
normahzed to Rplp2. Error bars represent the s.d. *, P<0.05. 
(TIFF) 



References 

1. Skinner MK (2011) Role of cpigcnctics in developmental biology and 
transgenerational inheritance. Birth Defects Res Part C Embryo Today Rev 
93: 51-55. doi;10.1002/bdrc. 20199. 

2. Cantone I, Fisher AG {201,3} Epigenetie programming and reprogramming 
during development. Nat Struct Mol Biol 20: 282-289. doi:10.1038/nsmb.2489. 

3. Bergemann AD, Cole F, Hirschhorn K (2005) The etiology of Wolf-Hirschhom 
syndrome. Trends Genet 21: 188-195. doi:16/j.tig.2005.01.008. 

4. Zollino M, Murdolo M, Marangi G, Pecile V, Galasso C, ct al. (2008) On the 
nosology and pathogenesis of Wolf— Hirschhorn syndrome: (jenot\'pe-phenoty]:)e 
correlation analysis of 80 patients and literature review. Am J Med (icnct G Se- 
min Med Genet 148G: 257-269. doi:10.1002/ajmg.c.30190. 

5. Nimura K, Ura K, Shiratori H, Ikawa M, Okabe M, et al. (2009) A histone H3 
lysine 36 trimethyltransferase links Nkx2-5 to Wolf-Hirschhom syndrome. 
Natiirc 460: 287-291. doi:10.1038/nahu-e08086. 

6. Andersen EF, Carey JC, Earl DL, Corzo D, Sutrie M, et al. (2013) Deletions 
involving genes WHSCI and LETMl may be necessary, but are not sufficient to 
cause Wolf-Hirschhorn syndrome. Eur J Hum Genet EJHG. doi:10.1038/ 
cjhg.201,S.192. 

7. Barron BJ, Robinson L, Tran HD, Lamki L {1994} Scintigraphic Manifestations 
of "Sternal Guppmg."J Nucl Med 35: 1 167-1 169. 

8. Mazzeu JF, Krepischi-Santos AG, Rosenberg G, Lourengo CM, et al. (2007) 
Widening the clinical spectrum of Rtt-Rogers-Danks/Wolf-Hirschhom syn- 
dromes. Genet Mol Biol 30: 339-342. doi:10.1590/S1415-47572007000300007. 

9. Komori T, Yagi H, Nomura S, Yamaguchi A, Sasaki K, et al. (1997) Targeted 
disruption of Cbfal results in a complete lack of bone formation owing to 
maturational arrest of osteoblasts. Cell 89: 755—764. 

10. Ghuang LSH, Ito K, Ito Y (2013) RUNX family: Regulation and diversification 
of roles through interacting proteins. IntJ Cancer 132: 1260-1271. 
doi:10.1002/ijc.27964. 

1 1 . Liakhovitskaia A, Lana-Elola E, Stamateris E, Rice DP, van't Hof RJ, et al. 
(2010) The essential requirement for Rurrxl in the development of the sternum. 
DevBiol 340: 539-546. doi:16/j.ydbio.2010.02.005. 

12. Kimura A, Inose H, Yano F, Fujita K, Ikeda T, et al. (2010) Runxl and Runx2 
cooperate during sternal morphogenesis. Development 137: 1159-1167. 
doi:10.1242/dcv.045005. 

13. Wang Y, Bclilower RM, Dong Y, Schwarz EM, O'Kcde RJ, et al. (2005) 
Runx 1 / Ami 1 /Cbfa2 Mediates Onset of Mesenchymal Cell Differentiation 
Toward Ghondrogenesis. J Bone Miner Res 20: 1624-1636. doi:10.1359/ 
JBMR.050516. 

14. Schroeder TM, Kahler RA, li X, Westendorf JJ (2004) Histone deacetylase 3 
interacts with rurrx2 to repress the osteocalcin promoter and regulate osteoblast 
differentiation. J Biol Chem 279: 41998-42007. doi:10.1074/jbc.M403702200. 

15. Lee M, Partridge NC (2010) Parathyroid hormone activation of matrix 
metalloproteinase- 1 3 transcription requires the histone acetyltransferase activity 
of p300 and PCAF and p300-depcndcnt acetylation of PCAF. J Biol Chem 285: 
38014-38022. doi:10.1074/jbc.Ml 10.142141. 

16. Kato S, Inouc K. Voun M-Y {2010} Emergence of the osteo-epigenome in bone 
biology. IBMS BoncKEy 7: 314-324. doiho. 1 1 38/20100464. 

17. Mikkelsen TS, Ku M, Jaffe DB, Issac B, Licbcrman E, et al. (2007) Genome- 
wide maps of chromatin state in pluripotcnt and lineage-committed cells. Nature 
448: 553-560. doi:10.1038/nature06008. 

18. Barski A, Guddapah S, Cui K, Roh T-Y, Schones DE, et al. (2007) High- 
Resolution Profiling of Histone Methylations in the Human Genome. Cell 129: 
823-837. doi:10.1016/j.cclL2007.05.009. 

19. Lucio-Eterovic AK, Singh MM, Gardner JE, Veerappan GS, Rice JC, et al. 
(2010) Role for the nuclear receptor-binding SET domain protein 1 (NSDl) 



Figure S2 Effects of electroporation of siRNA to 
MC3T3-E1 cells. Gene expression was examined by qRT- 
PCR and normalized to 18S. BMP3, 3 days culture with BMP2. 
BMP6, 6 days culture with BMP2. n = 3. Error bars indicate SD. 
*, p<0.05. N.S., not significant. 
(TIFF) 

Table SI Primers used in this research. 

(TIFF) 

Aclcnowledgments 

We thank M. Okado and T. Hiraoka for their technical assistance, K. Ura 
for providing mice, and the members of the GTS laboratory for their 
discussions and support. 

Author Contributions 

Conceived and designed the experiments: KN YK. Performed the 
experiments: YL KN WL KS. Analyzed the data: YL KN WL KS. 
Contributed to the writing of the manuscript: YL KN. 



methyltransferase in coordinating lysine 36 methylation at histone 3 with RNA 
polymerase 11 function. Proc Nati Acad Sci U S A 107: 16952-16957. 
doi: 1 0. 1073/pnas. 1 002653 1 07. 

20. Hu M, Sun X-J, Zhang Y-L, Kuang Y, Hu G-Q, ct al. (2010) Histone H3 lysine 
36 methyltransferase Hypb/Setd2 is required for embryonic vascular remod- 
eling. Proe Nati Acad Sci 107: 2956-2961. doi:10.1073/pnas.0915033107. 

21. Edmunds JW, Mahadevan LC, Gktyton AL (2007) Dynamic histone H3 
methylation during gene induction: HYPB/Setd2 mediates all H3K36 
trimedrylation. EMBO J 27: 406-420. doi:10.1038/sj.emboj. 7601967. 

22. Carrozza MJ, Li B, Rorcns L. Suganuma T, Swanson SK, et al. (200.5) Histone 
H3 Methylation b\ Srl2 Directs Dcacrtylation of Coding Regions by Rpd3S to 
Suppress Spurious Intragenic Transcription. Cell 123: 581-592. doi:10.1016/ 
j.ccll.2005. 10.023. 

23. Venkatcsh S, SmoUe M, Li H, Gogol MM, Samt M, et al. (2012) Set2 
methylation of histone H3 lysine 36 suppresses histone exchange on transcribed 
genes. Nature 489: 452-455. doi:10.1038/naturell326. 

24. Pinney DF, Emerson CP (1989) lOTl/2 cells: an in vitro model for molecukir 
genetic analysis of mesodermal determination and differentiation. Environ 
Healdr Perspect 80: 221-227. 

25. Czekanska EM, Stoddart MJ, Richards RG, Hayes JS (2012) In search of an 
osteoblast cell model for in vitro research. Eur Cell Mater 24: 1-17. 

26. Nimura K, Ishida C, Koriyama H, Hata K, Yamanaka S, et al. (2006) Dnmt3a2 
targets endogenous Dnmt3L to ES cell chromatin and induces regional DNA 
metirylation. Genes GeUs 11: 1225-1237. doi:10.1111/j.l365-2443.2006.010 
12.x. 

27. Takcichi .M, Nimura K, Mori M, Nakagami H, Kancda Y {2013} The 
transcription factors Tbxl8 and Wtl control the cpicardial epithclial- 
mescnchymal transition through bi-dircetional regulation of Slug in murine 
primary cpicardial cells. PloS One 8: c57829. doi:10.1371/journal.pone. 
0057829. 

28. Miyazono K, Maeda S, Imamura T (2004) Coordinate regulation of cell growth 
and differentiation by TGF-P superfamily and Runx proteins. Oncogene 23: 
4232-^237. doi: 10. 1038/sj.onc.l 207 131. 

29. Ducy P, Zhang R, Geoffroy V, Ridall AL, Karsenty G (1997) Osl2/Gbfal: A 
Transcriptional Activator of Osteoblast Differentiation. Cell 89: 747—754. 
doi:10.1016/S0092-8674{00)80257-3. 

30. Kitabayashi 1 (1998) Interaction and functional cooprralifin of the leukemia- 
associated factors AMLl and p300 in mvcloid cell differentiation. EMBO J 17: 
2994-3004. doi: 1 0. 1093/cmboj/ 17.11 .2994. 

31. Boumah CE, Lee M, Selvamurugan N, Shimizu E, Partridge NC (2009) Runx2 
recruits p300 to mediate parathyroid hormone's effects on histone acetylation 
and transcriptional activation of the matrix metalloproteinase- 13 gene. Mol 
Endocrinol Baltim Md 23: 1255-1263. doi:10.1210/me.2008-0217. 

32. Javed A, Gutierrez S, Montecino M, van Wijnen AJ, Stein JL, et al. (1999) 
Multiple Gbfa/AML Sites in the Rat Osteocalcin Promoter Are Required for 
Basal and Vitamin D-Responsive Transcription and Contribute to Chromatin 
Organization. Mol Cell Biol 19: 7491-7500. 

33. Yoh SM, Lucas JS, Jones KA (2008) The lwsl:Spt6:GTD complex controls 
cotranseriptional mRNA biosynthesis and HYPB/Setd2-mediatcd histone 
H3K36 methylation. Genes Dcv 22: 3422-3434. doi:10. 1 101 /gad. 1720008. 

34. Komori T (2006) Regulation of osteoblast differentiation by transcription 
factors. J Cell Biochem 99: 1233-1239. doi:10.1002/jcb.20958. 

35. Gnmstein M (1997) Histone acetylation in chromatin structure and transcrip- 
tion. Natiire 389: 349-352. doi:10.1038/38664. 



PLOS ONE I www.plosone.org 



10 



September 2014 | Volume 9 | Issue 9 | el 06661 



Whsci Fine-Tunes Bone-Related Genes 



36. Zhang S, Xiao Z, Luo J, Hr N, MahliosJ, ct al. (2009) Dosc-Dcpcndcnt Effects 37. Lou Y, Javcd A, Hus.sain S, Colby J, Frederick D, ct al. (2009) A Runx2 

of Runx2 on Bone Development. J Bone Miner Res 24; 1889—1904. threshold for the cleidocranial dysplasia phenotype. Hum Mol Genet 18:556— 

doi:10.1359/JBMR.090502. 568. doi:10.1093/hmg/ddn383. 



PLOS ONE I www.plosone.org 



11 



September 2014 | Volume 9 | Issue 9 | el 06661 



